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Abstract 

In this paper we present an 0(n log n) algorithm for finding a maximum flow in a directed 
planar graph, where the vertices are subject to capacity constraints, in addition to the arcs. 
If the source and the sink are on the same face, then our algorithm can be implemented in 
0(n) time. 

For general (not planar) graphs, vertex capacities do not make the problem more difficult, 
as there is a simple reduction that eliminates vertex capacities. However, this reduction does 
not preserve the planarity of the graph. The essence of our algorithm is a different reduction 
that does preserve the planarity, and can be implemented in linear time. For the special 
case of undirected planar graph, an algorithm with the same time complexity was recently 
claimed, but we show that it has a flaw. 

1 Introduction 

The problem of finding maximum flow in a graph is a well-studied problem with applications 
in many fields, see the book of Ahuja, Magnanti and Orlin [1] for a survey. The maximum 
flow problem is also interesting if we restrict it to planar graphs, which are graphs that have 
an embedding in the plane without crossing edges. The case of planar graphs appears in many 
applications of the problem, for example road traffic or VLSI design. The special structure of 
planar graphs allows us to get simpler and more efficient algorithms for the maximum flow and 
related problems. 

In the maximum flow problem, usually the arcs of the graph have capacities which limit the 
amount of flow that may go through each arc. We study a version of the problem in which the 
vertices of the graph also have capacities, which limit the amount of flow that may enter each 
vertex. This version appears for example when computing vertex disjoint paths in graphs, and 
in other problems where the vertices model objects which have a capacity. 

Ford and Fulkerson [31 Chapter 1. 11] studied this version of the problem. They suggested 
the following simple reduction to eliminate vertex capacities. We replace every vertex v with 
a finite capacity c by two vertices v' and v" . The arcs that were directed into v now enter v', 
and the arcs that were directed out of v now leave v" . We also add a new arc with capacity c 
from v 1 to v". Unfortunately, this reduction does not preserve the planarity of the graph [7]. 
Consider for example the complete graph of four vertices, where one of the vertices has finite 
capacity. This graph is planar. If we apply the construction of Ford and Fulkerson we get a 
graph whose underlying undirected graph is the compete graph with 5 vertices. This graph is 
not planar by Kuratowski's Theorem. 

The most efficient algorithm for maximum flow in directed planar graphs without vertex 
capacities, to date, was given by Borradaile and Klein [2] (Weihe |11] gave an algorithm with 
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the same time bound but assuming certain connectivity condition on the graph). Their paper 
also contains a survey of the history of the maximum flow problem on planar graphs. The time 
bound of the algorithm of [2] is 0(n log n) where n is the number of vertices in the input graph. 
Borradaile and Klein ask whether their algorithm can be generalized to the case where the flow 
is subject to vertex capacities. 

A planar graph is a st-planar graph if the source and the sink are on the same face. Hassin [I] 
gave an algorithm for the maximum flow problem in directed st-planar graphs without vertex 
capacities. The bottleneck of the algorithm is the computation of single-source shortest-path 
distances, which takes 0(n) time in a planar graph, using the algorithm of Henzinger et al. [5]. 

Khuller and Naor [7] were the first to study the problem of maximum flow with vertex 
capacities in planar graphs. They gave various results, including an 0(ny / Iogn) time algorithm 
for finding the value of the maximum flow in st-planar graphs (which can be improved to 0(n) 
time using the algorithm of [5]), an 0(n log n) time algorithm for finding the maximum flow in 
st-planar graphs, an 0(n log n) time algorithm for finding the value of the maximum flow in 
undirected planar graph, and an 0(n 1,5 log n) time algorithm for the same problem on directed 
planar graphs. If all vertices have unit capacities, then we get the vertex-disjoint paths problem. 
Ripphausen-Lipa et al. [TO] solved this problem in 0(n) time for undirected planar graph. 

Recently, Zhang, Liang and Chen [13] , used a construction similar to the one of [7] to obtain 
a maximum flow for undirected planar graph with vertex capacities that runs in 0(n log n) 
time. Their algorithm constructs a planar graph without vertex capacities, and then uses the 
algorithm of [2] to find a maximum flow on it, which is modified in 0(n log n) time to a flow in 
the original graph with vertex capacities. They also gave an 0(n) time algorithm for undirected 
st-planar graphs. Zhang et al. also ask in their paper if there is an algorithm that solves the 
problem for directed planar graph. 

In this paper we answer [2] and [13j . and show a linear time reduction of the problem of 
finding maximum flow in directed planar graphs with arc and vertex capacities, to the problem of 
finding maximum flow in directed graphs with only arc capacities. This problem is more general 
than the one for undirected planar graphs, since an undirected planar graph can be viewed as 
a special case of a directed planar graph, in which there are two opposite arcs between any pair 
of adjacent vertices. 

We show how to apply the constructions of [7J and [13] to directed planar graphs. Given 
directed planar graph, we construct another directed planar graph without vertex capacities, 
such that we can transform a maximum flow in the new graph back to a maximum flow in the 
original graph. Since the new graph does not have vertex capacities we can find a maximum 
flow in it using the algorithm of [2 J (or of [I] if it is an st-planar graph). The time bound of 
our reduction is linear in the size of the graph, therefore we show that vertex capacities do not 
increase the time complexity of the maximum flow problem also for planar graphs. 

In addition, we show that the algorithm of [13] unfortunately has a flaw. We give an 
undirected graph in which this algorithm does not find a correct maximum flow. Therefore, in 
fact our algorithm is also the first to solve the problem for undirected graphs. 

The outline of the paper is as follows: In the next section we give some background and 
terminology. In Section [3] we describe the construction of p] that we use, and in Section [5] we 
describe the one of [TO]. In Section [5] we characterize when a maximum flow in the constructed 
graph induces a maximum flow in the original graph. In Section [6] we show how to efficiently 
find such a flow in the constructed graph. Finally, in the last section we stick all the pieces 
together to get our algorithm. 
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2 Preliminaries 



We consider a simple directed planar graph G = (V,E), where V is the set of vertices and 
E is the set of arcs, with a given planar embedding. The planar embedding of the graph G 
is represented combinatorially, see [9] for survey on planar graphs. An arc e = (u, v ) G E is 
directed from u G V to v G V . We denote the number of vertices by n, since the graph is planar 
we have \E\ = O(n). 

A path P = (eo, e\, . . . , ek-i) is a sequence of arcs e^ = (ui,V{) such that for < i < k — 1 
we have Vi = Ui+\. If in addition v^—i = uq then P is a cycle. We say that a path P contains 
a vertex v, if either (u, v) or (v, u) is in P, for some vertex u. The path P = (eo, ei, . . . , e k -i) 
starts at no and ends at For t> G V, m(u) = {(u, f) | (u, u) G E} is the set of incoming 

arcs and out(v) = {(v, u) \ (v, u) G E} is the set of outgoing arcs. 

The graph G has two distinguished vertices, s G V is the source and t G V is the sirafc. 
The source s has no incoming arcs, and the sink t has no outgoing arcs. Every arc e£E, has 
a capacity c(e) > 0, and in addition every vertex v G V \ {s,t} has a capacity c(v) > 0. A 
capacity might be oo. We assume that the source and the sink have no capacities, if we wish to 
allow them to have capacities, we can add a vertex s' that will be the source instead of s, and 
an arc (s',s) with the desired capacity, and similarly add a new sink t' , and an arc (t, t') with 
the desired capacity. Note that this transformation keeps the graph planar, and even si-planar 
if it was so. It is easy to extend the given embedding to accommodate s' , t' , and the arcs (s' , s) 
and (t,t'). A graph without vertex capacities can be viewed as a special case in which c(v) = oo 
for every vertex. 

A function / : E — > R is a flow function if and only if it satisfies the following three 
constraints: 

< /(e) < c(e) Ve G E (1) 
/(e)<c(«) VveV\{s,t} (2) 

E /( e )= E /( e ) V^GFUM} (3) 

edin(v) edoutlv) 

Constraints ([I]) are the arc capacity constraints, Constraints j2|) are the vertex capacity con- 
straints and Constraints @ are the flow conservation constraints. 

We say that e G in(y) carries flow into v if /(e) > 0, and that e' G out{y) carries flow out 
of v if /(e') > 0. 

The uaZwe of a flow / is X^eemft) /( e )> the amount of flow which enters the sink. If the value 
of / is then / is a circulation. Our goal, in the maximum flow problem, is to find a flow 
function of maximum value. 

For a flow function / we define a cycle C to be a flow-cycle if /(e) > for every arc in C. 
We extend this definition to every function / : E — > R, even if it is not a flow. If a function / 
has no flow-cycles we say that / is acyclic. An acyclic flow is a flow function which is acyclic. 

Let e = (u, v) we denote rev(e) = (v,u). For a flow function /, we may assume that 
/ does not contain an arc e such that both /(e) > and f(rev{ej) > 0, because otherwise 
the flows in both directions can cancel each other. For a path P = (eo,ei, . . . ,ek-i) we let 
rev(P) = (rev(e k -i),rev(e k - 2 ), • • • , rev(e )). 

The planar embedding of G partitions the plane into connected regions called faces. For a 
simpler description of our algorithm, we fix an embedding of G such that t is on the boundary 
of the infinite face. It is easy to convert any given embedding to such an embedding [9]. 

The dual graph Gd of G has a vertex D(h) for every face h of G, and an arc D(e) for every 
arc e of G. The arc D (e) connects the two vertices corresponding to the faces incident to e. 
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Figure 1: A planar graph and its dual graph. The vertices of G are dots, and its arcs are solid. 
The vertices of Gd are circles, and its arcs are dashed. The bold arcs are an arc-cut in G and a 
cut-cycle in Gd- Capacities are not shown in this figure. 

The arc D(e) is directed from the vertex that corresponds to the face on the left side of e to 
the one of the face on the right side of e. Intuitively, Gd is obtained from G by turning the arcs 
clockwise. The dual graph Gd is planar, but it is may have loops or parallel arcs. Every face h 
of Gd corresponds to a vertex v in G, such that the arcs that bound h are dual to the arcs that 
are incident to v. See Figure HJ The capacity of e G E, c(e), is interpreted in Gd as the length 
of D{e). 

In the construction we present below we add undirected edges to directed graphs. Each such 
undirected edge uv can be represented by two antiparallel directed arcs (u,v) and (v,u), with 
the same capacity. If e is an undirected edge, then D(e) is also undirected. 

2.1 Residual cycles 

In this section we present the algorithm of Khuller, Naor and Klein [8] that finds a circulation 
without clockwise residual cycles, in a directed planar graph. The complete details of the 
algorithm can be found in [8] and also in [2]. We also present an extension of this algorithm 
that changes a given flow into another flow, with the same value, without clockwise residual 
cycles. We use this algorithm later to get the linear time bound for our reduction. This algorithm 
of [H] was given for directed planar graphs without vertex capacities, so for this section assume 
that the graph G does not have vertex capacities. 

In this section we also assume that if e S E then also rev(e) E E. This assumption can be 
satisfied, without changing the problem, by adding the arc rev(e) with capacity for every arc 
e such that rev(e) is not in E. 

Let / be a flow in G. The residual capacity of an arc e with respect to / is defined as 
c r (e) = c(e) — /(e) + f(rev(e)). In other words, the residual capacity of e is the amount of flow 
that we can add to e, or reduce from rev(e). The residual graph of G with respect to / has 
the same vertex set and arc set as G, and the capacity for each edge e is c r (e). A residual arc 
with respect to / is an arc e with a positive residual capacity. A residual path is a path made 
of residual arcs. A residual cycle is a cycle made of residual arcs. 

Given embedded planar graph G with arc capacities, Khuller et al. [H] showed how to find 
a circulation f r in G, such that there are no clockwise residual cycles with respect to f r . Their 
procedure is as follows. 

Let hoo be the the infinite face of G. Find the shortest path distances from D(hoo) to every 
vertex of Gd- Define a potential function <f> that assigns to every face h of G, the shortest 
distance from D(hoo) to D (h) in Gd- Let e be an arc of G with face hi to its left and face h r 
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to its right. If 4>{h r ) > <j){h() then f r (e) = 4>(h r ) — 4>(hi) (otherwise it is implied that f r (rev(e)) 
gets the value (f>(hi) — 4>{h r ) and f r (e) = 0). 

The function f r satisfies the two constraints of a flow [S] (without vertex capacities). The 
edge capacity constraints are satisfied, because for each e, the arc D(e) connects D(hi) to D(h r ), 
so <f>(h r ) ~ ^{ht) is at most the length of D(e), which is the capacity of e. The flow conversion 
constraints are satisfied because for every cycle in Gd, the sum of f r (e) — f r (rev(e)) for the arcs 
of the cycle is [HE], and the vertices which incident to a specific vertex in G form a cycle in 
Gd- Therefore, f r is a circulation. 

The correctness of this algorithm follows from the fact that any clockwise cycle C in G 
encloses some face h. The shortest path from D(hoo) to D(h) in Gd must contain an arc D(e) 
dual to an arc e of C. Because D{e) is in the shortest paths tree, the algorithm assigns to f r (e) 
the value of the length of D(e), which is the same as the capacity of e. Therefore the residual 
capacity of e with respect to f r is c(e) — / r (e) + f r (rev(e)) = 0, and e is not a residual arc with 
respect to f r . Therefore, C is not a residual cycle with respect to f r . 

The bottleneck of the algorithm of [8] is the computation of single-source shortest-path 
distances from the vertex Dfooo) in the dual graph. Henzinger et al. [5] showed how to find 
these distances in a planar graph in 0(n) time, so the algorithm of [8] can be implemented in 
the same time bound. 

Given a flow / in G, we wish to find a flow /' with the same value, such that /' does not 
have clockwise residual cycles. We can use the algorithm of [8j which we described above for 
this problem. 

Let G' be the residual graph of G with respect to /. We find a circulation f' r in G', such that 
G' does not have clockwise residual cycles with respect to f' r , using the algorithm of [5]. Define 
/' to be the sum of / and f' r , that is /'(e) = max{0, /(e) + f' r {e) — [f(rev(e)) + f' r {rev{e))]}. In 
other words, we add to /(e) the flow of f' r {e), and let the flows on e and rev(e) to cancel each 
other. 

The function /' satisfies the two constraints of a flow without vertex capacities. The capac- 
ity of an arc e in G' is c(e) — /(e) + f(rev(e)) and therefore f' r {e) is smaller than this capacity, 
therefore /(e) + f' r {e) — [f(rev(ej) + f' r {rev(e))} < c(e), so /'(e) < c(e) and the arc capacity con- 
straints are satisfied in /'. The conservation constraints are satisfied, because these constrains 
are satisfied for / and for f' r , and /' is the sum of these two flows. 

The value of the flow /' is the sum of the values of / and f' r . Since f' r is a circulation, its 
value is 0, and so the value of /' is the same as the value of /. 

The flow /' has the desired property that it has no clockwise residual cycles. To show that, 
we show that if C is a clockwise residual cycle in G with respect to /', then C is also a residual 
cycle in G' with respect to f' r , contrary to the way we find f' r . Let e be an arc of C, and 
assume for contradiction that e is not residual in G' with respect to f r . From our assumption 
f' r {e) = c r (e) = c(e) — /(e) + f(rev(e)) and f' r (rev{ej) = 0. Therefore, /'(e) = c(e) and e is not 
residual in G with respect to /', in contradiction to the fact that it is a member of C. 

Lemma 2.1. Let G be a directed planar graph without vertex capacities, and let f be a flow in 
G. We can find a flow /' in G, with the same value as f , such that /' does not have clockwise 
residual cycles, in 0(n) time. 

3 Minimum cut 

In a graph without vertex capacities, a cut S is a minimal subset of E such that every path from 
s to t contains an arc in S. To avoid ambiguity later, when we introduce cuts that may contain 
vertices, we call such a cut an arc-cut. See Figure[TJ The value of an arc-cut S is J2ees c ( e )- The 
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Figure 2: Construction of G c and G e for the graph in Figure [H The graph G c is presented as 
the dual graph of G e . The newly added (undirected) edges are without arrowheads. Capacities 
are not shown in this figure. 

minimum cut problem asks to find an arc-cut of minimum value. The fundamental connection 
between maximum flow and the minimum cut problems was given by Ford and Fulkerson [3] in 
the Max-Flow Min-Cut Theorem: 

Theorem 3.1. J3j/ The value of the maximum flow (in a graph without vertex capacities) is 
equal to the value of the minimum arc- cut in the same graph. 

Let C be a cycle in Gd- We say that C is a cut-cycle if it separates the faces corresponding 
to s and t, and goes counterclockwise around s (or equivalently, clockwise around t). See Figure 
[TJ The length of C is the sum of the lengths of its arcs. Johnson [6 J showed the following 
relation between the value of minimum arc-cut and the value of shortest cut-cycle: 

Lemma 3.2. Jfijl Let G be a directed planar graph without vertex capacities. Then the value of 
the minimum arc-cut of G, is the same as the length of the shortest cut-cycle in Gd- 

Ford and Fulkerson O Chapter 1. 11] extended the definition of cuts to graphs with vertex 
capacities. In such a graph, a cut S is a minimal subset of E U V such that every path from s 
to t contains an arc or a vertex in S. The value of a cut S is similarly defined as X^xeS c ( x )- 
Ford and Fulkerson also presented a version of the Max-Flow Min-Cut Theorem for graphs 
with vertex capacities, in this case the value of maximal flow (subject to both arc and vertex 
capacities) is equal to the value of the minimum cut (which contains both arcs and vertices). 

Khuller and Naor |7j extended Lemma f3. 21 using a supergraph G c of Gd which they construct 
as follows. Let h be a face of Gd that corresponds to a vertex v of G with finite capacity. We 
add a new vertex Vh inside h and connect it by an (undirected) edge of length c(v)/2 to every 
vertex on the boundary of h. See Figure [2J 

Lemma 3.3. VfjJ The values of the maximum flow and minimum cut in G are equal to the 
length of the shortest cut- cycle in G c . 

4 The extended graph 

Zhang, Liang and Jiang [12] and Zhang, Liang and Chen |13j construct the extended graph for 
an undirected graph with vertex capacities. We use the same construction for directed planar 
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graphs with vertex capacities. The extended graph is defined as follows. We replace every vertex 
v £ V which has a finite capacity with d vertices Vq, ■ ■ ■ , Vd-i, where d = \in{v)\ + \out(v)\ is the 
degree of v. We connect every to mod d with an (undirected) edge of capacity of c(v)/2. 

We make every arc that was adjacent to v, adjacent to some vertex V{ instead, such that each 
arc is connected to a different vertex «j, and the clockwise order of the arcs is preserved. We 
identify the new arc (u,Vi) or (vi,u) with the original arc (u,v) or (v,u). The resulting graph 
is denoted by G e , and the cycle that replaces v £ V in G e by C v . The graph G e is a simple 
directed planar graph without vertex capacities. The arc set of G e contains the arc set of G. 
See Figure [2j 

From the construction of G e and G c follows that G c is the dual of G e . Let v be a vertex 
with finite capacity and let h be the corresponding face in Gd- Then, in G e we replaced v with 
C v , and in G c we placed Vh inside h. The edges which connects Vh to the boundary of h are 
dual to the edges of C v . 

Combining Theorem 13.11 Lemma 13.21 and Lemma 13.31 we get that the value of the maximum 
flow in G e is the same as the value of minimum arc-cut in G e , which is the same as the value 
of the shortest cut-cycle in G c , which equals to the value of the maximum flow in G. And the 
next lemma follows. 

Lemma 4.1. The value of the maximum flow of G is equal to the value of the maximum flow 
ofG e . 

5 Reduction from the extended graph to the original graph 

Let / be a flow in G e , we define / to be the restriction of / to the arcs of G. The next lemma 
generalizes and corrects the result of Zhang et al. \13\ Theorem 3]. 

Lemma 5.1. Let f be a flow function in G e . If f is acyclic then f is a flow function in G. 

Proof. We show that / satisfies all three conditions that a flow function in G should satisfy. 

Since the capacities of common arcs of G and G e are the same, / clearly satisfies arc capacity 
constraints. 

Let v G V. If c(v) = oo then v is also a vertex of G e , with the same incident arcs, and 
therefore / satisfies the flow conservation constraint at v. Otherwise, the amount of flow that 
enters C v in /, is the same amount that enters v in /. Also, the amount of flow that leaves 
C v in / is the same amount that leaves v in /. Therefore we obtain that / satisfies the flow 
conservation constraint at v by summing up the flow conservation constraints that / satisfies 
for the vertices in C v . 

It is left to show that / satisfies the vertex capacities constraints. (Note that these con- 
straints are irrelevant for / since in G e we do not have vertex capacities.) Let v £ V be a 
vertex with finite capacity. First, we show that the acyclicity of / implies that in the cyclic 
order around v of arcs with positive flow that are incident to v in G, the arcs carrying flow into 
v are consecutive, and the arcs carrying flow out of v are consecutive. Note that this claimed 
consecutiveness is restricted to arcs with positive flow, so arcs e with /(e) = can appear 
anywhere in the cyclic order of the arcs incident to v. 

If there is only one arc carrying flow into v then the claim is trivial. Otherwise, consider two 
arcs e and e' carrying flow into v. We show that if we cyclically traverse the arcs incident to v 
clockwise starting from e, we either traverse all arcs carrying flow out of v before traversing e', 
or we traverse them all following e' but before we get back to e. Since this hold for every pair 
of arcs e and e' carrying flow into v, the desired consecutiveness follows. 
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Figure 3: The path R from u to t must cross one of the paths Q or Q' from u to v, and thus 
creates a flow-cycle. The shaded area is H. 



Since e carries flow, there is a path P from s to v of arcs with positive flow that ends with 
e. Similarly, there is a path P' of arcs with positive flow from s to v that ends with e'. 

Let u be the last vertex of P, before v, that also appears on P'. The vertex u must also be 
the last vertex on P', before v, that also appears on P, as otherwise we get that / contains a 
flow-cycle. Let Q be the suffix of P that starts at the arc of P that goes out of u, and let Q' be 
the suffix of P' that starts at the arc of P' that goes out of u. (These arcs are uniquely defined 
since / does not contain a flow-cycle.) Since both Q and Q' goes from it to v the arcs of Q and 
Q' partition the plane into two regions, denote them by H and H'. 

The vertex v is obviously not t since v has a finite capacity. Furthermore, since there is an 
arc with positive flow outgoing of each vertex along P and P', none of these vertices can be t. 
Therefore t is either inside H or inside H'. We assume without loss of generality that t is in H'. 

Let d = (v, w) be an arc that carries flow out of v, such that w £ H. There must be a path 
R that starts with d and carries flow from v to t. See Figure [3j Since w is in H and t is in H', 
we get that there exists a flow-cycle that starts with a prefix of R and ends with a suffix of Q 
or of Q' , contradicting the assumption that / is acyclic. Therefore, all arcs that carry flow out 
of v, carry it to a vertex in H'. Our claim that the arcs carrying flow into v are consecutive 
and the arcs carrying flow out of v are consecutive then follows. 

An arcs e with /(e) > incident to v in G corresponds to the same arc e with /(e) > 
incident to a vertex in the cycle C v in G e . Therefore among the arcs incident to C v with /(e) > 
those that carry flow into vertices in C v are consecutive. Therefore, we can find two edges b 
and b' of C v such that if we remove them C v splits into two parts, such that all arcs that carry 
flow into vertices in C v are incident to one of these parts, and all arcs that carry flow out of 
vertices in C v are incident to the other part. See Figure HI All the flow which enters C v must 
go through b and 6', in order to leave C v , because of the flow conservation constraints on the 
vertices of C v . The total capacity of b and b' is c(v), and therefore the total flow that enters C v 
in / is at most c(y). Thus it follows that the total flow that enters v in / is also at most c(v), 
and the vertex capacity constraint at v holds. □ 



6 Canceling flow-cycles 

In order to use Lemma 15. II we must find a maximum flow / in G e such that / is acyclic. In this 
section we show how to do that. 

The algorithm of Zhang et al. |13|, Section 3] for undirected planar graphs finds a flow / in 
G e and than cancels flow-cycles in / in an arbitrary order. They call the resulting flow f a and 
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Figure 4: The edges b and b' separate between the incoming flow into C v and the outgoing flow 
from this cycle. Dotted arcs do not carry flow. The arrowheads of edges of C v indicate direction 
of flow. 




s t 



Figure 5: A counterexample to the algorithm of [13]. The edges of the original undirected graph 
G are solid. The vertex v has capacity 1, the edges of C v are dotted. The flow in every solid edge 
is 1, the flow in every dotted edges is 1/2, in the specified direction (the edges are undirected). 
The bold edges form a flow-cycle in G, after we cancel it we remain with an acyclic flow in G, 
but the amount of flow that enters v is 2. The correct solution is to cancel the flow in the two 
internal flow-cycles. 

claim that this flow satisfies vertex capacities constraints. This approach is flawed. Figure [5] 
shows an example on which the algorithm of [13] fails. After we cancel flow-cycles in / in an 
arbitrary order it is possible that there is no flow /' in G e whose restriction to G is /' = f a , 
and therefore Lemma 15. II does not apply. 

As the example in Figure [S] shows, it is not enough to cancel arbitrary flow-cycles in /. 
We can cancel a flow-cycle in / only if there is a cycle C in G e that contains it, such that we 
can reduce flow along the cycle C. In this case the cycle rev{C) in G e is a residual cycle with 
respect to /. Therefore, in order to cancel a flow-cycle in G with respect to / we must cancel a 
residual cycle in G e with respect to /. Canceling a arbitrary residual cycle is not enough, since 
we always want to reduce the flow that / assigns to arcs, and never to increase it. 

Let / be a flow in G e . We define a new capacity function d on the arcs of G e which 
guarantees that the flow in arcs of G never increases beyond /. For e 6 E we let d(e) = /(e). 
The arcs of G e which are not in G are arcs of C v for some vertex v, for these arc we do not have 
to limit the flow to the amount in /, so we set d(e) = c(e) = c(v)/2. The flow function / is also 
a flow function in G e with the new capacity d, by the way we defined d . Since c'(e) < c(e) for 
every arc e, every flow in G e with capacity d is also a flow in G e with the original capacity c. 

Instead of canceling the residual cycles one by one, we apply to G e and d the algorithm in 
Section [2.11 and find a new flow /' with the same value as /, such that there are no clockwise 
residual cycles in G e with respect to /' and d. The following lemma shows the crucial property 
of/'. 

Lemma 6.1. The restriction f of f to G does not contain counterclockwise flow-cycles. 

Proof. Assume, for a contradiction, that there is a counterclockwise flow-cycle C with respect to 
/' in G. We choose C such that C does not contain any other counterclockwise flow-cycle inside 
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its embedding in the plane. We show that we can extend rev(C) to a clockwise residual cycle 
with respect to /' in the graph G e with capacity c', in contradiction to the way we constructed 
/'• 

For every arc e of C, /'(e) > 0, and therefore rev{e) is a residual arc with respect to /' and 
d . If C does not contain a vertex v £ V with c(v) ^ oo then rev(C) is a clockwise residual 
cycle with respect to /' and c', and we obtain a contradiction. 

Let v be a vertex in C with c(v) 7^ 00. Let (u,v) and (v,u') be the arcs of C which are 
incident to v. These arcs correspond to arcs (u,Vi) and (vj,u') in G e , where Vi and Vj are in 
C v . Let P be the path from Vj to Vi which goes counterclockwise around C v (recall that C v is 
undirected in G e ). To show a complete residual cycle in G e , we argue that P is a residual path 
in G e with respect to /' and c', so we can use it to fill the gap between Vj and Vi in rev(C). An 
arc e of P is not residual if and only if /'(e) = c'(e) = c(v)/2. Without loss of generality we 
assume that the vertex Vk in the path P is followed by Vk+i- 

Let ej = (uj-i, Uj) be the last arc in the path P from vj to Vj. Assume for contradiction 
that &i is not residual with respect to /' and c'. Then /'(e^) = c(v)/2 and so the total flow 
which enters into Vi in /' is Yle&nOvi) f'( v i) — f'((. u i v i)) + f'( e i) > c ( v )/^- The only remaining 
arc that can carry flow out of is (^,1^+1). Because /' satisfies flow conservation constraints 
f'((vi, Ui+i)) > c(v)/2. But this is impossible since the capacity of the arc (i>i,Ui+i) is c(v)/2. 
Therefore, /'(e^) < c{v)/2 and e^ is residual with respect to /' and c' . 

We now proceed by induction. Assume by induction that we already know that the arc 
e fc+i = (vk,Vk+i) on the path P from Vj to «j is residual with respect to /' and c'. If k = j 
then we are done. Otherwise, we prove that e^ = (^fc-i,^) G P is also residual with respect 
to /' and c'. Since ek+i is residual it follows that /'(efc+i) < c(v)/2 . Let e' be the single arc 
of -E 1 incident to v^. If e' is directed out of and /'(e') > then since C is a cycle and e' 
is inside C in the embedding of G, there must be a path carrying flow that starts with e' and 
continues to another vertex on C. (Recall that t is incident to the outer face.) This implies 
that there is a counterclockwise flow-cycle with respect to /' and G inside the embedding of C, 
in contradiction to the choice of C. Therefore e' does not carry flow of /' out of Vf. in G e . This 
implies, by the conservation constraint on Vk, that /'(e^) < f'(ek+i) < c(v)/2, so e^ is indeed 
residual with respect to /' and c'. 

We showed that there is a residual path from Vj to V{. Since v was an arbitrary vertex with 
c(v) > on C it follows that we can extend rev(C) to a residual cycle in G e with respect to /' 
and c'. Since C is a counterclockwise cycle, the residual cycle we got from rev(C) is a clockwise 
cycle. This contradicts the definition of /', and therefore a counterclockwise flow-cycle C with 
respect to /' and G does not exist. □ 

We repeat the previous procedure symmetrically, by defining a new capacity c" which re- 
stricts the flow in G to the flow in /', and applying a symmetric version of the algorithm of 
Section \2. 11 This way we get from /' a flow /" of the same value, such that /" does not contain 
clockwise flow-cycles in G. For every e G E we changed the flow such that /"(e) < /'(e) < /(e), 
so we did not create any new flow-cycles. Therefore we have the following lemma. 

Lemma 6.2. The flow function /" has the same value as the flow function f . The restriction 
/" °f f" to G is acyclic. 

7 The algorithm 

Combining together the results of the previous sections we get an algorithm for finding maximum 
flow in a directed planar graph with vertex capacities. 
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First, we construct G e from G by replacing each vertex that has a finite capacity with C v as 
defined in Section HI Next, we find a maximum flow / in G e , which is a directed planar graph 
without vertex capacities. Last, we change / to another flow /" as in Section 

According to Lemma f6.2l the flow /" is maximum flow in G e , and its restriction /" is acyclic. 
By Lemma [5. 11 the function /" is a flow in G. And by Lemma [4.11 /" is a maximum flow, since 
the amount of flow that /" carries into t is the same as the amount of flow that / carries into t. 

The construction of G e from G takes 0(n) time. The computation of /" from / also takes 
0(n) time using the algorithm we described in Section f2. II Therefore, the only bottleneck of our 
algorithm is finding /, a maximum flow in a directed planar graph without vertex capacities. 

Theorem 7.1. The maximum flow in a directed planar graph with both arc capacities and vertex 
capacities can be computed within the same time bound as the maximum flow in a directed planar 
graph with arc capacities only. 

The algorithm of Borradaile and Klein [2] finds a maximum flow in directed planar graph 
with arcs capacities in 0(n log n) time. If G is a st-planar graph, then G e preserves this property. 
In this case the algorithm of Hassin [1] , using the algorithm of [5| for single-source shortest-path 
distances, finds a maximum flow in 0(n) time. 
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